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Приводится методика расчёта радиальных подшипников сколь- ТБе 1есблиаие оЁ са1сайие фе гад1а]! зПтае Беагиез оЁ 


жения бесконечной и конечной длины, работающих на страти- шйпие апа Апие еп гаппше оп а доц е-1ауег згаййе4 


фицированном двухелойном смазочном материале. Методика бнсаиЕ 15 ргезещед. Те 1есбоате 15 Базе оп Фе изе оЁ 


построена на использовании автомодельной переменной, что фе зе!Ё-зипа]ае уама е Ша аПо\уз оМашише Ше ехаси 


позволяет получить точное автомодельное решение задачи как в р о 
р 7 р а . а - зе!-зппПаг рго ет зо[айоп Бо ш ро[аг ап ш суппанса1 
полярных. так и в цилиндрических координатах, а также значе- ый. 
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рованных слоёв. Приведённые задачи усложнены одновремен- р | р 
и — _ и й Бу Фе зппаЦапеоч$ соп$14егайоп оЁ фе абисапЕ у15созКу- 
ным учётом зависимости вязкости смазочного материала от дав- 
ления. наличием адаптированного к условиям трения опорного ргеззиге тейайоп, Ше ргезепсе оЁ ше Елейоп-аЧаре4 зир- 


профиля подшипниковой втулки и влиянием осевой подачи сма- рогит8 рго@е оР Фе Беаппз РизВ, ап4 Ше еНесе оР Ве ах1а| 


зочного материала в подшипнике конечной длины. Получены бисапе Еее4 ш Фе Нойе-еп2 Беаппо. Аз а гези, сот- 


расчетные модели Для двухслойных расслаивающихся смазоч- 
ных материалов, численный анализ которых позволил уетано- 
вить влияние переменных факторов на основные рабочие харак- 
теристики подшипников — компоненты вектора несущей спо- 
собности, силы трения и оптимальные значения параметра 
опорного профиля, подачи смазочного материала, 6его вязкостно- 
го соотношения в стратифицированных СЛОЯХ. 


Ключевые слова: Радиальный подшипник, двухслойный сма- 
зочный материал, осевая подача смазочного материала, зависи- 
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функции тока. автомодельное решение. 


Введение. Как известно [1-3], при наличии в смазочной жидкости частиц присадок или продуктов износа, а также за 
счет пристенной адсорбции и ориентации молекул, ВОЛиИЗи опорной поверхности ПОДШИПНИКОВОЙ втулки происходит 





стратификация смазочного материала на слои с различной вязкостью. Слоистое течение вязкой несжимаемой жидкости в 
зазоре упорного и радиального подшипников рассматривалось в работах [4-9]. Существенный недостаток существую- 
щих методик заключается в том, что в расчетной модели не учитывается зависимость вязкости от давления. При боль- 
птих значениях давления в смазочном слое вязкость смазки существенно возрастает и возникает необходимость учета 
зависимости вязкости от давления. 

Постановка задачи. Цель работы — разработать аналитический метод расчета радиальных подшипников, работающих 


БИр://уезби К. доп ва.га 


на двухслойной смазке, с учетом зависимости ее вязкости от давления. Также необходимо оценить влияние вязкостного 
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= ниннианый 
отношения слоев; параметра, характеризующего границу раздела слоев на основные рабочие характеристики подшипни- 
ка; параметра, обусловленного зависимостью вязкости от давления при осевой подаче смазки. 

Задача № 1. Рассматривается течение двухслойного смазочного материала в зазоре радиального подшипника беско- 
нечной длины. Вал вращается с угловой скоростью 62, а подшипниковая втулка с адаптированным профилем опорной 
поверхности неподвижна. Зависимость вязкости от давления выражается формулой 


; 


ре”, (1.1) 
В полярной системе координат с полюсом в центре вала (рис. 1. 1) уравнения контуров вала, границы раздела 
слоев и адаптированного контура опорной поверхности запишутся в виде: 
с: =К: а:И=и +694 + 69050 — ©4$шт 50: 
со: г’ =, +есоз0-— Ат 00, где @аЕ[0,1|. (1.2) 


Здесь го — радиус вала; 7› — радиус подшипника; го + ба — радиус границы раздела слоев; параметры А и © 
характеризуют адаптированный профиль опорной поверхности. 





Рис. 1.1. Схема радиального подшипника с двухслойной стратификацией 
жидкого смазочного материала 


Исходные уравнения и граничные условия. В качестве системы ИСХОДНЫХ уравнений берется безразмерное уравне- 


НИе ДВИЖЕНИЯ вязкой несжимаемой жидкости с учетом зависимости вязкости от давления и уравнение неразрывности 
г 
@е ы ми е РА. 
[6 Я9 


где размерные величины 7’, и,, ®,, р’, и, связаны с безразмерными и, и;,о;, р, и; следующими соотношениями: 


ар. ов ВИ Во р (1.3) 
ОГ 09 








„! — ” я / — ” / — я / — / — р — ” ” 

= о’, 5=ОО, и = и, р=р.р, ШЕНож, 9=-ц, (1.4) 
где д“ — экспериментальная постоянная; Но; —_ характерные вязкости смазочных слоев, и, ? о; — компоненты вектора 
скорости. © — угловая скорость вращения вала, р, — характерное давление. 


Система уравнений (1.3) решается при следующих граничных условиях: 














и ‚=7(0) — 0, 12> „=7(0) — 0: и ‚=ай— и ГЕМЙ> я г=ай — 92| „=ал> 
р} 2} 
во) Мо? СО И; 
1 МИНИ 2 2 1 Е ! 
я, [7=ай = ПА, |^=ай? | зал =0" ©). 
@), Но 9 ; 
й(9) =1+1с0$0 —-п 3100, п=е/б, п! =/А/б. (1.5) 





Точное автомодельное решение. Точное автомодельное решение системы уравнений Сы удовлетворяющее гра- 
НИЧНЫМ УСЛОВИЯМ (1.5% будем искать. ИСПОЛЬЗУЯ функцию Тока, предполагая, что поле скоростей и давлений в сма- 
очном СЛое является потенциальным: 


САИ, 


и; = 
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— би (,,. а. 
т (7(г.О). ОЕ А. + Г. (г.О). АИ; — р; (5). 
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А. й 5 а с 
о о. (1.6) 


е® 4 17(0) 130)’ 


ы И > = = М> = 
ее. 
Мо Мо 


Осуществляя подстановку (1.6) в (1.3) и (1.5). получим: 


— — 
— 


и =С, 9=а, ш+е=0, =6, 5=а, +65 =0, (1.7) 
$ (0)=0; #(0)=0, 5 (0)=Ё ч,()=0, #050, 


5>(1) —- # 1 (06) -- > (0%), 5; (00) № 5> (0%). 1 (0%) = И (04). 





[оз ] 
26, Ми), [5+ [55 ©МЕ-0 (1.8) 
НЯ НЯ 0 © 


Решение задачи (1.7)-(1.8) находится с помощью непосредственного интегрирования. В результате получим: 


ЕО ее 
{= 6, -— +5. +6., 9 =СР-+СЕ+С 
2 2 
2 2 
Е - ы и. - ь 
АР = СС. РС: 125 == +555 Сб 
й 2 
=З =2 ЕЗ =2 
д = 6+0, №2 =-4 2-26-01, 
3 2 3 


Ане = лей 26.1.6456] лье = Алей 9.6) +8,- 6] 


0 
к (9) = | 


0 


49 
(+ с0$09-пу $т ‹90)^_ 


(1.9) 


м м 


Для определения постоянных с,(1=2,3,...,11) с, с, 4, с› используется следующая алгебраическая система 
из 14 уравнений с 14 неизвестными: 
Ст —= 1. Со — 0. СЗ — 0. 


> | — — 
— — 


| | | 
-<--се-+е1=0 с-чеч+о =0 сЪ-+е+с.: =0 
1 $ 1 › “] $169 > 62] 4165 : 
3 2 2 2 








выд Ц > ее Ц > вы о (2%) 
С — = С 6 С> — = С> ь С> о ти ь 
НЯ НЯ /3(2т) 
= — мо ыы = У 
И НЯ 
2 2 
— Я тая 0 
2 2 
2 2 


— 


— 9 т 
аси -<-—-с0-—со =0, 
Г" т, 9 =%. 


3 2 3 2 
гы 06 9? 6? © о | 1 
С —— +6 — +619 -<— < — -бо“+а-+с-—+со =0. (1.10) 
6 Р. 6 6 


С учетом (1 4) гидродинамическое Давление р определяется следующим выражением: 


р=1+] 1+5 | Зло ИЛИ р=1+] 1+5 | -лзво (1.11) 
ОЛА < РА 


Записывая систему (1.10) в матричной форме, получим: 
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М -х=Ь, (1.12) 

где х=Е: СД; С55 С с р = 0,0, — 60: 0;—2} 

ие. 2 2 0 0 
з(2п) 
1 0 0 2 2 
М =| -о?-+1 0 0 З/о-3а“`+3 6-6а. 

2 /2(2%) 

(1 К)о =2 е--Т: = 0 0 
Л3(2п) | 

“(К - 1) 0 0 2—1) 2 


Решая систему (1.12), получим 





ь м т) (3 — 607 — Во? +30" — о В" +320“ — бро”) 
= _ б-+6Ао” — 607 ме (т) 
д ры (ак — м + ПД 
ОТ За? — За +303 +3 ба? — З% +3932 — За? 4? — браз +304) 
ре 1. (2п) 
^ (А -м- ПА 
4—40? + 4ког — 4 + 4Ко? — ЗКог” + Зо” +1 
= б= о СНА, СЕКС, 
д д 
А =-407 +10" — бКо” + 4Ког + К^ое* + 4Коь — Ко” — 40% + бо, 
О ЕТ | 
ыы Пи С То, 
/3(2п) 2ло В 2 
я #1, 
пы (1.13) 
5 (2т) 


Основные рабочие характеристики подшипника. Нормализованные показатели расходов От И (9 двухслойной 
смазочной жидкости определяются выражениями 
ог о” ее (97 о” 
6 2 6 2 6 2 


Используя формулы (1.13) и (1.14) для компонент нормализованной поддерживающей силы и нормализован- 


— 


НОЙ СИЛЫ трения, придем к следующим выражениям: 





а р. и Г: 050 — а ен лы | с0$(® — 12-1 А со$ (6 + 1/2 —1 | (97 | 
ГоРе 








0 ы: Л! 20 Ф+1 Ф-1 2) 
Г К, Ге: 1020 ст: | зш@®-102л эзщ(о®-+12л | ь. о 
„= й И = С — — Пезыиа- би: довел оков ПНА: БН — 
` РэТо ь 20 АЛ. 20 ф-—1 +1 2) 
Тс 2 
— ГЙ % ТИ 1% ы 
[р = —- а ге”? Ю (1.15) 


| 126) 100 
ОР 0 Й (0) И ) ё=0 


Итоги численного анализа полученных аналитических выражений для базовых эксплуатационных характери- 
стик узла трения показывают следующее: 

— при значении вязкостного параметра © Е [0, Пнесущая способность достигает максимума при величине па- 
раметра контура профиля ® = 5; 

— с увеличением вязкостного параметра а несущая способность подшипника возрастает на 5—10 %. 
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0.3 ы 


Рис. 1.2. Зависимость вертикальной компоненты безразмерной несущей способности Ау от параметра эксцентриситета 1 и 
параметра адаптивного профиля <: [-а=0.1: 2-а= 0.2; 3 —а= 0,5; Ч-а= 0,9 


Задача №2. Рассматривается установившееся течение двухслойной смазки в зазоре радиального подшипника конечной 
ДЛИНЫ. Вал вращается с ПОСТОЯННОЙ УГЛОВОЙ скоростью (2, а подшииниковая втулка неподвижна. Подача смазочного 
материала осуществляется в осевом направлении, а его вязкость выражается следу ющей экспоненциальной зависимо- 
Стью ОТ Давления: 


: 


Шо”. (2 


ы Подшипникавая втулка „ „ 


Шейка вала 





Рис. 2.1. Расчетная схема 


В цилиндрической системе координат (у’, 0, 2’)с началом в центре вала уравнения контуров вала, границы 


раздела слоев и адаптированного контура опорной поверхности запишутся в виде. 


со: =: < :Г = +60 + 926050 — 9Аз100: 
Я (2.2) 
с.г =, +е605$0- Азто0. где чеЕ[0,| = И. 


Исходные уравнения и граничные условия. Движение смазочной среды описывается системой уравнений ДвВиИжЖе- 
ния вязкой несжимаемой жидкости для приближения ТИПа «ТОНКИЙ слой» И уравнениями неразрывности 


22 





- -,б их „9 би. 00, ди 
_ 5 - ле О\ = ле В АО (1=1,2), (2.3) 
бр д С 0 9 6% 
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где размерные величины р’, р. и, о, р; и; ‚ описывающие распределение скоростей и давлений в смазочном слое. 


связаны с соотвегствующими стандартизированными величинами и, 2, и;, ©;, р;, И; соотношениями 


и РЕ и т —: ПИ й Г и __ и __ 
м=з +0, г =и2, ч=Оцм, х=Фцо, ии =Оби,, р=р,р, М.=Нойу. (2.4) 
Здесь и;, 5, м — компоненты вектора скорости смазочных сред; р’ — гидродинамическое давление в сма- 
зочных слоях. и, — динамические коэффициенты вязкости в смазочных слоях; @ — экспериментальная постоянная; 
57 
ОР: 
А =——, р. — характерное давление. 
2? 28 
Нос 270 


Граничные условия на поверхности вала и подшипника записываются в виде 
и; |4 0=0. © |-0=1 20.2) = р0л.2); "|. 0=0, р(0.0)=1 ро = Рк 
Рз 
Ио |“ =05 02|)» М2 |--= 0. (2.5) 
На границе раздела слоев граничные условия записываются в виде: 
и би = И [590% 1 р ой = №2 |. ой ТТ ве и | 50а 
в а То 
де ЕСИ" и Ор гЕСЙ> бя г=ой— и 6» 
1(0) =1+1с0$0 —\, $100. ц=е/д, п = 4/6. (2.6) 


Граничные Условия (2.5) означают прилипание смазки к поверхности вала и подшипника, периодичность гид- 


И: 
. и. ее ой (09), 
о; 


родинамического Давления, а Также то, что при 9=0 давление задано в двух сечениях. 

Условия (2.6) означают равенство скоростей, касательных и Нормальных напряжений на границе раздела слоев. а 
Также условие существования слоисТтого течения смазки, Т. ©. требу ется, чтобы скорость Точек границы раздела слоев в 
каждой точке была направлена по касательной к контуру раздела слоев. 





Точное автомодельное решение. Точное автомодельное решение системы уравнений (2.3), удовлетворяющее гра- 
ничным условиям (2.5) и (2.6), ищется в виде 














_ О; Ор Е и 
И +©,(г.9)> ©;= +И(г.9 Ау; = 4, (5). м; =". (&,0), 
56 (7,0), 5, 5, (7,9). и; = У) (5,9) 
ы о р 
(/(г,0)=-и, (2) (0), У(г.0=5,). &=—. 
| 7 | | >” ^ 9 
| = о | р а - 
№ е % = 80| _ Е - -_ 0-+(—4)а2, Ле “? = ( 8] я — Ю +(—0)а52, 
„И (0) (0) ” то тЫ И г(0) | 
12 и. | к 
а = а> Е. а> = (р -), а= 2 а. (2.7) 


ы р. Нот 
Подставляя (2.7) в (2.3), (2.5) и (2.6), придем к следующей системе дифференциальных уравнений и гранич- 
ных условий к ним: 
И = >, 5 =9, т +50 = 0, \> =>, 5> =. 


5 


— 








> = ай” (9), 2 — ай” (0), и НЕ. = 0; (2.8) 
0 Е 


$ (0)=0: #(0)=0, З®=Ь #(0)=0, $, @=0. 
и-(1)=0, 55(10=0. м>(10=0. 5 (9) =5,(9). Ш) =и, (9), 











ве _№я = _ ИЯ Но2 
51 (©) =—= 55 (4%). (0) = 2 (о), а=—а, ©0=-“с,, а=—ба,, 
Ни Я Я Я Нот 
2) а {9 1 
> Ом; ь О» лы И" 
и (©, 0) = +. (м,0), — а = го и а | 5 (©) + | 5 (2) =0. (2.9) 
ОЕ. м 58 . 
у (е 


Учитывая, что разделение смазочного материала на слои происходит вблизи неподвижной твердой поверхности, 
т.е. при значениях ©, близких к единице, условие раздельного течения смазки (ий,(00)/5,(00)) = ой (9) в принятом нами при- 


ближении удовлетворяется. На самом деле из граничного условия 
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1 
1 (0) + 05 (©) + | 5. (ЕМЕ=0 (2.10) 
[9 
следует 
| 
— У ЕЛЕ 
— Н> (0) (64 
(а ча + | 227 |= 0. (2.11) 
02(0.) 02(0.) 
[92 
Используя теорему о среднем значении, будем иметь 
Е > 
№ г -| 5 (0) 1, «9 \ ия : * , | 
0. (0) 0, =0, а = (0,1). (2.12) 
т, 0; о, 0, 
Так как 5.(@ )<5.(9), (1-@)<<1, следовательно, с точностью до членов 0| —2^И-— с) |, будем иметь 
_ 0510) 
9. 1 
| 5, (4 =0, | 5. (24 =0. 
О © 
Решение задачи (2.8)-(2.9) находится непосредственным интегрированием. В результате будем иметь 
22 | 22 
, и В Е = рее. р = ры = 
\И = <> —-+<6+с;. в =а--+с+с, м =ай’--часч+а., 
2. 2 2 
2 — г —” _ #2 
из #2 дед 9 #2 
> _ хз > м > > > 


р =8.1, (0) +2. (0)+а2+Ь, р =6.1,(0)+6,./.(@)+а2+Ь,, 


0 


Л. (0) = | НИЕ. ИИ —. (2.13) 
> (1+1 с0$09 — 1 т <0) 


Для определения постоянных с;({=2,...,11), аА(у=\2,3,4), а,6,а,с> решим следующую алгебраическую систему 


из 18 уравнений с 18 неизвестными: 


Ст — 1. Со — 0. СЗ — 0. — С г — С ть. - С 1 аи 0. С а - С - Со — 0. 
3 2 
у . | 
> 9 = 09 = 1 
с. —— +620 +6 =6———=640-65 =0; © +6; =0, 
| й. й 
7 5 — 
к 0 о. м о = С а2(2т) 
++ -<---6щ-в =0 З=АС, © =. О-о, 
] 6 Я $ о ] 1 2 2 | 
(т) 
3 2 3 2 
к В (0 В 0% (0 +5 1 й И о 
а + чея-а--&-ф-топ+а-+о-+о = 
| 6 7 | 8— о 1 с о : 
6 2 6 2 
+; = (а, +с:), ба + с = (+ са), 
ь ь 


— — 


— 


4. =0, ат ++, =0, ай? — +4а = ай? —- +4. +а4.. 


ао +а, = Е (а>й^а +4.), КЕМ./ Ш. (2.14) 
вы 


Решение системы (2. 14) сводится к решению следующих двух матричных уравнений: 


М-Х=Ь. (2.15) 


Бр: //уезши К. Чопа .га 
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Ре 2 2 0 0 
/3(2т) 
| 0 0 2 2 
М= ю-в+1 0 0 Зо’ -3а`+3 6-6“ 
о У>(2п) 
але ^^ О > 0 0 
92(к—1) 0 0 20-1 2 
№-У=т, (2.16) 
| о? р” 
где У= Ч; Ч; 44, т = | - — (а -а); Пока фа); фа р 
и -( -| 
МЕТ -А 0. (2.17) 
от 1 


Решая матричные уравнения (2.15) и (2.16), с учетом (2.14), будем иметь 








ь - № (3— бо == о + ЗЕ^о — бо”) 
= _ б6+6ко” — 60/7 Ри 30) 
| А `“ (© -—и+ ПА 
= а — 39° — 3% +3+ бк” = ЗК +За К” — За” К” — бког + Зо) 
РЕ От) 
^ (иА-&-+ПА | 
4—4? +40? — 40 + За” +4Ког — Зко” +1 (2.18) 
с = бб =, СЕКС, СЕКС, 
А А 
А=-—407 +1+0” — бо” +4Ког + К^о +4 — Ко” — 4 +69”, 
г ь 2 | ре рые | 
ое = Ч (соз 2п°-П. с = 51 + и (с0$2ло® — , 
/3(27) 2% - 2% 
с Ека. р = =. р. 4. =0, 4 =-а.-а ай 
- От) ‚ | я. 
1 А-а | 
— Ка>о — а 0 + а>(1-=& 
5—> Ша;а - па+а, 1—0) 


4 = | 4, = Кай“ + ЧК — ай “о. 
Ко — 9-1 | 


Определение ОСНОВНЫХ рабочих характеристик подшипника. Для гидродинамического Давления в принятом нами 
приближении приходим к следующему уравнению. 
29 2012 

— р — У 1 __ 


0. 
А: А 





ор —= 2 р 2—0 


„т 


б = м = Го — 
Решая это уравнение, с точностью до членов 0(1.0.). $ п будем иметь 


1 и) а? С. 
=1+—|<../.(9)+55.7. (0) 1+— [+2 1+— |. 2.19 
р А, [< 2(0)- соЛз( ]| > ие > (2.19) 


Приведем гидродинамические силы, возникающие в смазочных слоях и воздействующие на вал, к центру вала 
(7. Для ортогональных проекций безразмерного главного вектора этих сил, действующих на единицу длины вала, 
главного момента и момента сопротивления осевому движению, Получим выражения 
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Выражения для расходов смазочной жидкости в окружном и в осевом направлениях. записанные в Нормали- 


зованной форме, будут иметь следующий вид: 


— 


—- 





Рис. 2.2. Зависимость величины безразмерной несущей способности 


—— 


Ол 





2п 1 


[Ио + [504 | 


0 [02 


ъ 2 
В ФУ й 


) : 


(2.20) 


7 от параметра расслоения смазочного материала п 


и параметра адаптированного опорного профиля в при значениях вязкостных параметров: / = а = 0,1, А/=1; 2-а= 0,3, К = 


1; 3 —а-0.5, А2/К1 = 1,2; 4- а = 0.9. А/К = 1,3 


Проведение качественного анализа полученных аналитических выражений (2.20) для основных эксплуатаци- 


онных характеристик узла трения позволяет сделать следующие выводы: 


— компонента те безразмерной поддерживающей силы существенно зависит от вязкостного параметра а и 


вязкостных отношений А>/К1; 


— В рамках принятых предположений максимум несущей способности достигается при значении параметра 


1 
опорного профиля ® = я : 


— учет зависимости вязкости от давления приводит к повышению расчетной несущей способности. 


— В рамках принятых предположений, с увеличением значений вязкостного параметра а, сила сопротивления 


для осевого движения смазочного материала увеличивается. 


Рассмотрев расчетные модели радиальных подшипников можно сделать следующие выводы: 


- анализ полученных расчегных моделей радиальных подшипников © круговым опорным профилем, работа- 


ющих на вязком несжимаемом смазочном материале при его двухслойной стратификации, позволил установить, что 


как обычные Подшипники, так и Подшипники © Пористым покрытием на шейке вала имеют зависимость несущей спо- 


собности от параметра вязкости смазочного материала 4, вязкостного соотношения его стратифицированных слоев К2 


БИр://уезби К. доп ва.га 
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и параметра а. разделяющего слои. Причем с ростом этих параметров величина несущей способности подшипника 
увеличивается, а при значении & = 0,2 имеет максимум: 

- расчетные модели конечных и бесконечноразмерных радиальных подшипников, с адаптированным к усло- 
виям трения опорным профилем и двухслойным стратифицированным вязким смазочным материалом при зависимо- 
сти его вязкости от гидродинамического давления, показывают при анализе, что несущая способность подшипников 
увеличивается с ростом вязкостного параметра а, параметра | и параметра опорного профиля ®, при значении кото- 
рого ® = 0,5 несущая способность имеет максимум; 

- Теоретические исследования модели радиального подшипника, смазываемого материалом с двумя страти- 
фицированными слоями, один из которых является вязким, с зависимостью вязкости от гидродинамического давле- 
ния, а другой обладает вязкопластичностью, позволили установить, что несущая способность подшипника растет с 
увеличением параметра вязкости @, параметра стратифицированных слоев @ (при а = 0,2 имеет место максимум) и 
параметра пластичности А. 
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